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A model for a Heisenberg antiferromagnet on a pyrochlore lattice with exchange and dipole-dipole 
interactions is studied via mean-field theory. In treating the dipoles by use of the Ewald method, 
a soft (critical) mode with a unique ordering wave vector is selected for all strengths of the dipole- 
dipole coupling. For weak dipoles a partially ordered, three sublattice spin structure (P state) , with 
Qord = ^5 51 is selected. A fully ordered, four sublattice spin structure (F state), with q — 000, 
competes with the P state and becomes the stable structure as the temperature is reduced. Our 
results are compared against other theoretical calculations and connection to recent experimental 
results for Gd2Ti207 are discussed. 
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Frustration is ubiquitous in condensed matter and 
arises whenever a system cannot minimize its total clas- 
sical ground state energy by minimizing the energy of 
its pairwise interactions, pair by pair. In the context 
of magnetism, frustration commonly occurs in triangu- 
lar or tetrahedral unit based lattice structures with an- 
tiferromagnetic (AFM) interactions, geometric frustra- 
tion, and prevents the magnetic moments (spins) from 
pointing antiparallel to each other. Two and three di- 
mensional lattices of corner sharing triangles or tetrahe- 
dra with nearest neighbor AFM Heisenberg exchange are 
highly frustrated and are particularly interesting. There, 
mean-field theory (MFT) finds a massive degeneracy of 
soft (critical) modes with no unique preferred state cho- 
sen at some critical temperature As a consequence, 
frustrated systems exhibit a propensity for thermal and 
quantum fluctuations 0, 0| ■ The possibility that exotic 
novel quantum states arise in frustrated quantum spin 
systems is currently attracting much attention. 

It is generally expected that perturbations to the near- 
est neighbor spin Hamiltonian select a unique ordered 
state. In the insulating pyrochlore systems R2M2O7, 
where is a magnetic trivalent rare earth ion (R=Gd, 
Tb, Ho, Dy), and M=Ti^+ or Sn^+, non-magnetic 
(Fig. ^), the leading perturbation to the nearest- 
neighbor Hamiltonian is the long-range magnetic dipole- 
dipole interaction. Due to its (i) long-range 1 /r^ and (ii) 
anisotropic r ■ S spin— space coupling nature, one would 
naively expect dipolar interactions to select a unique and 
robust soft magnetization density (mode) at an ordering 
wave vector Qord- However, explicit calculations in two 
specific cases have found dramatically different behavior. 

In the Ho2(Ti,Sn)207 and Dy2Ti207 pyrochlores with 
local (111) Ising spins, it has been found that the 
dipolar interactions are to a very large extend self- 
screened 0- As a result, the system does not develop 
long-range order at a critical temperature Tc Dnn, 
where Dnn ^ 2 K is the dipolar energy scale at nearest- 



neighbors. Instead, the system develops an extensive 
quasi-degenerate "spin ice" state at T ^ D^n analogous 
to the proton disordered state in common ice water with 
the same residual extensive Pauling entropv0,0|. Inter- 
estingly, detailed calculations (mean-field |2] and Monte 
Carlo III) find that the dipolar interactions do not pro- 
duce an exact symmetry and thus allows for the weak 
selection of an ordered state at ordering wave vector 
Qord = 001. Monte Carlo simulations employing nonlocal 
loop dynamics find a first order transition, removing all 
residual entropy, to the Qord = 001 long-range ordered 
state at the critical temperature Tc/Z)„n ~ 0.07 ^ 1 V]. 

In the Gd2(Ti,Sn)207 pyrochlores, the Gd'^^ ions 
have S — 7/2 and L — 0, hence negligible single-ion 
anisotropy. They are, therefore, good realizations of an- 
tiferromagnetically coupled Heisenberg spins with weak 
long-range magnetic dipole-dipole interactions (20% of 
nearest- neighbor AFM exchange) 8]. Here too, the ef- 
fect of dipolar interactions are far from obvious. Two 
independent mean-field (Gaussian approximation) calcu- 
lations have found that long-range dipole-dipole inter- 
actions truncated at large distance (i.e., 100 nearest- 
neighbors) give rise to a line of degenerate soft modes 
along the cubic (111) direction Palmer and Chalker 

argued that including quartic terms in the theory leads 
to the selection of a four sublattice Neel ordered state 
at Qord = 000, see Fig. QJd. Experiments provide a dif- 
ferent picture: (i) Neutron diffraction measurements at 
T = 50mK find a partially ordered phase at q = 
see Fig. with one disordered sublattice [lO(. (ii) 
Mossbauer experiments |lll | find strong evidence for a 
paramagnetic (PM) fluctuating spin down to 35mK, too. 
(iii) Speciflc heat measurements 12J find two transitions 
at T = 0.7K and T — l.OK. All these experimental results 
are seemingly incompatible with the prediction of Ref. |^ 
In view of these perplexing results from theory and ex- 
periment, we have revisited the mean-field calculations 
of Refs. HandU We have found that the degeneracy of 
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soft modes along the q = 111 direction reported inl^and 
1^ is only approximate. Instead, we find that long-range 
dipolar interactions, when properly re-summed to infinite 
distance, give rise to a weakly selected unique soft mode 
at Qord — with one disordered sublattice, similar to 
what is found in experiments (Fig. ^) 10]. Just 0.7% 
below that first transition, we find a second transition to 
the q = 000 state of Fig. ^p. We also find that weak 
second and third nearest neighbor exchange only acts to 
suppress or enhance the transition temperature of the 
q = 000 mode relative to that of the qord 
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state. 
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FIG. 1: (a) The pyrochlore lattice as a fee Bravais lattice 
with a four atom basis, (b) The fully ordered q = 000 equal 
moment spin structure (F state) about a single tetrahedron. 
Each spin is parallel to an opposite side of a tetrahedron, each 
tetrahedron has zero net moment, ^iLi = 0- ("=) The 



partially ordered q ■ 



spin structure (P state) projected 



onto the plane perpendicular to the [111] direction. The solid 
(dashed) lines represent a tetrahedron projected out of (in 
to) the plane of the paper. The three sublattices with equal 
moments lie in a plane (kagome) and sum to zero with the 
fourth interstitial site, encircled, paramagnetic. There is a 
reversal of spin direction between spins in adjacent layers. 



We consider a model for large S-spins on a py- 
rochlore lattice with exchange and dipolar interac- 
tions but negligible single-ion anisotropy. The ap- 
propriate Heisenberg Hamiltonian, H, is given by 

H = -iE,,Ea,6E„,„J.t(^j)5r'"^f, where in- 
dices denote fee Bravais lattice points, (a, 6) 

describe the tetrahedral basis, and {u, v) define the 
components of classical Heisenberg spins, S'f'". The 
Juvihj) interaction matrix contains exchange and 
dipolar contributions and is given by Juvihj) — 



JnS; 



D. 



dd 
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nearest neighbor distance, Ddd = ^^nn with D = 
/j,o/x^/47ri?^,j, ^ is the magnetic moment occupying the 
pyrochlore lattice site, i?nn is the nearest neighbor dis- 
tance, and R'^j is the vector separation between spins 
and Sj. The nearest neighbor exchange, Ji, sets the 
energy scale (Ji < is AFM). The effects of weak long- 
range exchange (J2 and J3) are also considered. 

We first investigate the soft-mode spectrum of H by 
studying the momentum dependence of the static sus- 
ceptibility, x{q) wfthin MFT [HQ- The mean-field 
free energy is = Tr{piJ} + TTr{plnp} H 0, where 
Tr is a trace over spin configurations and p is the mean- 
field (one-particle) density matrix. The expression for T 
is expanded to quadratic order in mf, the vector order 
parameter, Fourier transformed, and then diagonalized 
by a normal mode transformation. The diagonal form of 
!F is used to calculate the correlation function {Sf ■ Sj) 
in the definition of x(<7). 
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where Jjv^ij, Afim restricts the exchange to the iV-th 



where P — 1/T in units of fee = 1 and A'ccii is the num- 
ber of Bravais lattice points, and q is wave vector in the 
first Brillouin zone. The final form for the q-dependent 
susceptibility is x('Z)/^ccii = iP/5)Ec.,^\F^iQ)\^i^ ~ 
X'^{q)/3T)~^ , where the Greek indices label the normal 
modes (a — 1,2,3,4 and ^ — 1,2,3) at each value of 

q, and F^iq) = Ea,uUi:^i<l) = EaU^'"iQ) is a 3- 
component vector. U (q) is the unitary matrix that di- 
agonalizes J^{q) in the sublattice space with eigenval- 
ues X'^{q). The static susceptibility is a maximum at 
the ordering wave vector, qord, which also defines the 
maximum eigenvalue, i.e., A™'''^(qord) = maxq{A™'*'^(q)}, 
where X™^'^{q) is the maximum eigenvalue at q and maxq 
selects the global maximum for all q. The mean-field crit- 
ical temperature is given by, = A™'^''(qord)/3. 

To find qord for our model, we search for the maxi- 
mum in the soft mode spectrum (A™'''^(q)). We consider 
AFM first neighbor exchange, Ji = —1.0, and a vari- 
able dipole-dipole interaction strength, D. The dipolar 
interactions in J7(q) are summed directly to a cut off 
separation distance, Rc, measured in units of the nearest 
neighbor distance and to the infinite range limit through 
the use of the Ewald method ^3] ■ The regime of weak 
dipole-dipole interactions is of greatest interest to us be- 
cause the value D/\Ji\ = 0.20 faithfully approximates 
Gd2Ti2O70. At this point, we find a spectrum in the 
(hhl) plane that appears degenerate along the {hhh) di- 
agonal, in agreement with earlier works (Fig. 6 in 
Ref. 8). However, when we examine the spectrum along 
the {hhh) diagonal in greater detail, we observe the selec- 
tion of a global maximum as one increases the cut off dis- 
tance of the dipolar lattice sum. Figure El displays results 
for several values of Rc and the Ewald limit. An order- 
ing wave vector occurs at qord = \ W (in units of 27r/a), 
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and emerges from a smoothly varying spectrum along 
(hhh) in the Ewald limit. We note that for dipole-dipole 
interactions summed to Rc < 500, X™'^^{q) displays sub- 
stantial variations making it difficult to ascertain the ex- 
istence of a critical ordering wave vector. These fluc- 
tuations are damped with increasing Rc and approach 
the Ewald limit for Rc > 1000. A similar sensitivity of 
the soft-mode spectrum to the dipole-dipole cut off has 
been observed in the spin ice materials In a broader 
expanse of phase space, we find three regions for an or- 
dering wave vector: (i) qord = for D/\Ji\ < 5.60, 
(ii) qr„,d = 001 for 5.60 < L»/|Ji| < 5.85, (iii) ^ord = 000 
for D/\Ji\ > 5.85. Earlier work that imposed a cut off 
on the dipole-dipole lattice sum found a degenerate spec- 
trum along the (hhh) diagonal for D/\Ji\ < 5.70, with 
the selection of an ordering wave vector near (000) for 
larger values of D, (see Fig. 2 in Ref. 1^. 

The spin structure for <7ord — h^^h can be obtained 
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from the corresponding eigenvectors, U^'"{^^^). About 
a single tetrahedron one infers that three sublattices have 
equal moments with all spins confined to a plane and 
a fourth sublattice that is PM, i.e., zero moment, see 
Fig. When this configuration is extended over the 
pyrochlore lattice and modulated by qord, one has or- 
dered kagome planes along the [111] direction with a re- 
versal of spin direction in alternating planes, all intersti- 
tial sites are PM. This partially ordered three sublattice 
structure (P state) has been observed in neutron scatter- 
ing experiments 01 and is supported by Mossbauer[lH 
experiments on Gd2Ti207. The mean-field transition 

\ max /111 
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temperature to this state is 
0.94008|Ji| iD/\Ji\ = 0.2). 

The mode q — 000 corresponds to a fully ordered 
(four sublattice) equal moment structure (F state) ini- 
tially predicted to be the ground state of a D/\Ji\ = 0.20 
model 1^ . We find the limit of stability for the paramag- 
netic state against an ordering to the F state occurs at 
Tj = A'°'''^(000)/3 = 0.93923 1 Ji I, which is very close to 
rP = 0.940081 Ji I 15]. Given the close proximity of 
and TjF at the Gaussian level, we study the evolution of 
the respective free energies as a function of temperature 
below both values of T^. 

The real-space mean-field theory (RSMFT) begins 
with the mean-field decomposition (5*°'" — m^'")(5'^'" — 

m^'") — applied to H. The free energy is determined 
from / = -(1//3) InZuF, where Zmf is the RSMFT par- 
tition function. For our model, one obtains 

1 



a,u b.v 



i,j a.b u,v 



i.a 



InZf, (2) 



where = (47r//3|Bf |) sinh(/3|Bf |) is the single-site 
partition function and is the local field. The sub- 
lattice magnetization m° (order parameter) is obtained 
from differentiation of / with respect to local field. 
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One finds m" = 
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where the 



FIG. 2: Soft-mode spectrum along (hhh) as a function of 
dipolar cut off, Rc, with _D/|Ji| = 0.20. (a) The spectrum 
develops Rc dependent local maxima for Rc > 1. (b) The 
fluctuations in X'^'^'^{hhh) are smoothed out as Rc increases 
with the Rc ^ <x limit obtained with the Ewald method. In 
the Ewald limit, the critical mode falls at Qord = 



field at site {i,a) is due to all the other moments in 
the system acting at site (i,a), i.e., Bf = J2u-^'i '" = 
J2u J2j.b,v ^™ It is convenient to factor out the 

magnitude of the order parameter and write, = mbf, 
where mf = mfif and fif is a unit vector that defines 
the orientation of the order parameter at site (i,a). Us- 
ing this form for Bf, one expands the right-hand-side of 
the above expression for m° for small m to obtain Tc for 
RSMFT, T^^ = |6^|/3. 

We work on a lattice of L x L x L cubic cells, or 16L^ 
spins, with periodic boundary conditions and D/\Ji\ — 
0.20. The lattice is tiled with one of the above order 
parameters. We consider systems of size L = 4, 6 and 
calculate the dipole-dipole interactions via the Ewald 
method in real-space |l4|. Our procedure is to solve 
self-consistently for the sublattice magnetization, m°, at 
temperature T/\Ji \ and then calculate the corresponding 
free energy from Eq.[21 i.e., f /\Ji\. The real-space transi- 
tion temperatures are = \bf^''\/3 = 0.94002|Ji| and 
= \bf^^\/3 = 0.939171 Jil for L = 6, in good agree- 
ment with the above q-space results. For L = 8, the 
real-space and q-space results for and Tj agree within 
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10"''. Free energy versus temperature curves for the F 
and P ordered states on a L = 4 lattice are displayed 
in Fig. 131 At high temperatures, near either or Tj, 
the free energy curves appear degenerate. However, the 
insets to Fig.|31demonstrate that the respective curves ac- 
tually cross. We performed linear fits to points from the 
two curves about the crossing temperature. Solving these 
equations we find the temperature Tcross ~ 0.93358 |Ji|. 
Above Tcross, the P state has a lower free energy and is 
thus the preferred state. Below Tcross, the free energy of 
the F state drops below that of the P state and eventu- 
ally becomes the stable ground state of the system. We 
find T = energies per spin of = — 1.40856 |Ji| and 
= -1.05737|Ji|. Hence, in this very simple mean- 
field picture, the three sublattice P state is the initially 
preferred structure at high temperatures, but the nearby 
four sublattice F state competes with it and very rapidly 
wins out as the temperature is lowered. 




T/|J,| 

FIG. 3: Free energy of the P and F order parameters for a 
L = 4 lattice with D/|Ji| = 0.20. At high temperatures the 
free energy curves appear to lie on top of one another with 
the F state clearly preferred at lower temperatures. The insets 
show the free energy above (P mode favored) and below (F 
mode favored) a crossing temperature Tcross/|Ji| ~ 0.93358. 

The proximity of the P and F states in MFT suggests 
it would be difficult to realize their experimental differ- 
entiation (or in Monte Carlo simulations). However, the 
inclusion of weak long-range exchange in our model 
separates the transition temperatures of these two states 
and effectively reduces their competition. We note that 
a J2 or J3 has no effect on the transition temperature of 
the P state ■ ^^e Fourier transform of the J2 and J3 
matrices, the particular combination of the lattice sym- 
metry and the q = tttttt phase yields zero matrix elements 
for J2 (tttttt) and Jy,{'K'K'K) in the above calculation of xio) 
at q = TTTTTT. Hence, we find that J3 = 0.02Ji (AFM) 
suppresses the F state to 0.859237| Ji| (T = 4). A 

small FM J2 has a similar effect but of reduced magni- 
tude. Conversely, a small AFM J2 or FM J3 enhances 



the F state by raising its transition temperature above 

c • 

Hence, experimentally reasonable values for long-range 
exchange |l6| could allow for the detection of two well- 
separated transitions in Gd2Ti207. Although two tran- 
sitions have been observed in specific heat measurements 
at T IK and T « 0.7K neutron scattering results 
only find the three sublattice P state down to T w 50mK 
|l8| . We note that a structure is purposed experi- 
mentally for the P state, but a 4-fc spin state is a possi- 
bility. In our RSMFT, a A-k spin structure corresponds 
to a unit cubic cell in which one of the tetrahedra is fully 
disordered If is unclear at this time whether it is 

a material (sample) specific effect or something else that 
causes the P state to freeze out in the neutron experi- 
ments as the material is cooled from T > 1 K down to 
50mK . Recent ESR results have also found peculiar 
behavior in Gd2Ti207, where a large [111] anisotropy 
develops spontaneously below 80K |20j. Clearly, more 
experimental and theoretical work is needed to elucidate 
the peculiar thermodynamic properties of Gd2Ti207. 

In conclusion, our calculations above have identified 
another system, apart from spin ice, where the long-range 
1/r^ nature of dipolar interactions cause a large degree of 
symmetry restoration on the q— dependent susceptibility. 
However, this symmetry restoration is incomplete, and 
gives rise ultimately to a well defined ordering state upon 
cooling from the fully disordered paramagnetic state. 
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